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0.1. INTRODUCCION 

Einstein nos presentó hace ya casi 100 años una descripción de la gravedad 
como geometría. Ahora, el Cálculo de Regge nos presenta una aproximación 
a la Relatividad General en la que la variedad Espacio-Tiempo es modelada 
mediante bloques 4-dimcnsionales con geometría plana ( Lorentziana ) llamados 
Simplex. Más exactamente, un complex simplicial geométrico cuadridimensional 
es utilizado como un modelo del espacio-tiempo. 

Este conjunto de notas es una pequeña introducción al método desarrollado 
por Regge^l y su utilización en la solución de problemas en Relatividad General. 

En este modelo, las propiedades geométricas de la variedad son expresadas 
en términos de la longitud de las aristas del complex. De esta manera, el Cálculo 
de Regge se convierte en una teoría discretizada en la que las longitudes de las 
aristas son las variables dinámicas fundamentales. 

Supongamos ahora que conocemos la geometría de una hipersuperficie como- 
de-espacio de la variedad 4-dimcnsional espacio-tiempo. Ya que las ecuaciones 
geometrodinámicas de Einstein son deterministas, la geometría a partir de este 
instante estará completamente predeterminada. 

En términos del Cálculo de Regge, esto quiere decir que sera posible calcular 
la longitud de las aristas de los nuevos simplex a partir de las longitudes de 
los simplex anteriores. Una vez se halla discretizado o esqueletizado toda la 
variedad, tendremos un "catálogo" con la longitud de todas las aristas o huesos, 
y con ello conoceremos la geometría completa de cada hipersuperficie. 



Parte I 

Cálculo de Reg 



Capítulo 1 



Simplex, Complex 
Simpliciales y Poliedros 



Comenazaremos definiendo los bloques fundamentales con los que se mode- 
lará una variedad n-dimensional. 

Deflnicin 1 Dados n+1 puntos enl m , (m > n), enumerados por po,Pi,P2, 
si estos puntos son geométricamente independientes (i.e no existe un hiperplano 
(n-1) -dimensional que contenga a los n+1 puntos) , entonces definimos el sim- 
plex n-dimensional T n =< pop\....p n > como : 



Así, tenemos que un 0-simplex Tq =< po > es un punto o vértice, un 1- 
simplex Ti = < P1P2 > es una linea o arista, un 2-simplex T% =< P0P1P2 > es 
un triangulo con su interior incluido y un 3-simplex T 3 —< P0P1P2P3 > es un 
tetrahedro sólido. 

Nota 1 Ya que T n es un subconjunto cerrado de R m entonces T n es compacto. 

Deflnicin 2 Sea T n =< poP\...p n > un n-simplex y < r < n un entero. Si 
tomamos r + 1 puntos Pi Q ,Pi t ^ ■■■tPí^. del conjunto Po,Pi, ■■■,Pn entonces estos 
forman un r-simplex T r =< pi a pi 1 ...pi r > y se denomina una r-cara del simplex 




donde c¿ > 



y 




(1.1) 



T 



Nota 2 El numero de r- caras de un n- Simplex es 
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Definicin 3 Sea K un conjunto finito de simplex en W 1 . Si estos simplex están 
unidos entre si y satisfacen: 

1. '&) 1) Si a es una cara de un simplex en K, entonces a € K 

2) Si T y T son simplex de K, entonces T C\T' es vacio o una cara 
deT y de T 

entonces K se denomina un Complex Simplicial. 

Nota 3 La dimensión de un Complex K se define como la máxima dimensión 
de los simplex del conjunto K. 

Definicin 4 Sea K un Complex. Si cada simplex perteneciente a K es consi- 
derado como un subconjunto de R™ (n > dim K) entonces la unión de todos los 
simplex es un subconjunto de IR". Este subconjunto se denomina el Poliedro \K\ 
del complex simplicial K. 

Nota 4 La dimensión de \K\ como subconjunto de W 1 es la misma de K. 
{áim\K\=á\mK). 



Definicin 5 Sea X un espacio Topológico. Si existe un Complex K y un ho- 
meomorfismo f : \K\ — ► X, entonces se dice que X es triangulable y el par 
(K,f) se denomina una Triangulación de X. 



1.1. Poliedros y variedades bi-dimensionales 

Comenzaremos ahora con una descripción inicial de la teoría de curvatura 
intrínseca en poliedros en un caso simple tal como lo es una variedad bidimen- 
sional (superficie) para luego realizar su generalización a n dimensiones. En 
general una superficie arbitraria puede ser considerada como el límite de una 
sucesión de poliedros cuyo número de caras va incrementándose. La aproxima- 
ción que nos ofrecen estos poliedros es mala si se observa de cerca la superficie, 
pero para un observador lejano esta es satisfactoria. 

1.1.1. Curvatura Integral Gaussiana 

En cualquier superficie puede definirse una curvatura integral Gaussiana 
mediante el uso de triángulos construidos con geodésicas. 

Supongamos un triangulo t cuyos ángulos internos son a, (3, 7. En gene- 
ral, para geometrías no-Euclidcanas tenemos a + (3 + 7 7^ n. Así, definimos la 
Curvatura Integral Gaussiana eat del triangulo t como: 



e Gt =a + [3 + "f-ir 



(1.2) 
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De esta manera, por ejemplo, en una esfera, la curvatura integral Gaussiana 
será: ect = j¡? donde A t es el área del triangulo t y R es el radio de la esfera. 

Ahora bien, si t se reduce hasta convertirse un punto P y el límite lím ^-p- = 

K(P) existe y es independiente de la forma de realizarse, se toma la función 
K(P) como la definición de Curvatura Gaussiana Local en el punto P. 

En el caso de la esfera tenemos K(P) = . 

De esta manera tenemos que la Curvatura Integral Gaussiana se puede es- 
cribir como: 

e Gt = J K{P)dA (1.3) 

T 

Sin embargo debemos hacer claridad en que la curvatura Integral puede 
existir bajo pocas suposiciones, mientras que la curvatura local puede no existir 
como una función ordinaria, sino como una medida. 

Si consideramos ahora un poliedro M, observamos que si tomamos un trian- 
gulo t ubicado completamente en una de las caras del poliedro, tendremos 
£g* = 0. Este mismo resultado se obtiene si el triangulo no incluye ninguno 
de los vértices de M. De aqui concluimos que K{P) = si P no es un vértice. 

Por otro lado, si t incluye el vértice V de M, la curvatura integral no de- 
penderá de la forma de t. Es decir que ect = £v solo depende del vértice que 
se toma. Para encontrar la deficiencia del vértice V tomamos la suma de todos 
los ángulos internos de las caras de M con vértice V, y esta debe valer 2ir — e v . 

En general si t contiene varios vértices V, W, X,etc la curvatura será ect — 
£v + £w + £x + ••■ y por lo tanto K(P) resulta ser una función Delta de Dirac 
con los vértices como soporte. 



Formula De Euler para Poliedros 

Para observar un ejemplo del tratamiento de la curvatura, consideremos 
un poliedro M compacto (i.c. finito y cerrado). El teorema para la curvatura 
integral de Gauss-Bonnet se puede escribir como: 

^Te„ = 2^(2 - N) (1.4) 

n 

donde N es la clase de M { e.g. N = para la esfera, N — 2 para el toro). 
Ahora bien, si suponemos, sin pérdida de generalidad, que las caras de M son 
triángulos, y denotamos por cr/„ el ángulo interno de la cara / con vértice n, 
tenemos por la definición de la sección anterior: 



f 



(1.5) 
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donde la sumatoria se realiza sobre las caras / que tengan a n como vértice. 
Tenemos además que 

n 

donde la sumatoria se realiza ahora sobre los vértices de la cara /. 

Ahora bien, la sumatoria doble ~^2<Jf n puede ser obtenida de dos maneras 

fn 

diferentes: Sumando primero sobre los vértices y luego sobre las caras tenemos: 

fn f n f 

donde / es el número de caras de M. 

Por otro lado, sumando primero sobre las caras y utilizando el Teorema de 
Guass-Bonnet tenemos: 

fn n f n n 

donde V es el número de Vértices de M. Ya que estas dos dos maneras de 
realizar la sumatoria deben llevar a resultados iguales tenemos: 

f = 2(V-2 + N) (1.9) 

Ya que las caras de M son triángulos tenemos que 2e = 3/ donde e es 
el número de aristas de M. Reemplazando este valor en la ecuación anterior 
obtenemos la formula de Euler para poliedros: 

V-e + f = 2-N (1.10) 

1.2. Matriz de Conexión y Transporte Paralelo 

Veamos ahora como una variedad M puede ser modelada mediante una red 
simplicial o complex. Una vez construida esta red, la geometría intrínseca de M 
vendrá dada por la longitud de sus aristas y por la llamada Matriz de Conexión, 
la cual es en escencia una lista de todas las caras, aristas y vértices de M junto 
con su relación mutua. Esta matriz nos permitirá conocer cuales vértices y 
aristas pertenecen a una cara dada, y por lo tanto en ella esta implícita toda 
la información topológica de M. Además, al conocer la longitud de todas las 
aristas, tendremos también las deficiencias e n de cada vértice. 

De esta manera, la noción de red simplicial reemplaza la red de coordenadas 
que se utilizan en la variedad, mientras que el tensor métrico será reemplazado 
por la longitud de las aristas. 
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Cuando el número de vértices y/o caras de M se incrementa, la Curvatura 
Local se puede aproximar a una función continua: K(P) = ge , donde g es la 
densidad de vértices y e su deficiencia. Esta aproximación es valida en el caso 
en que este producto varié suavemente dentro del triángulo t. 

Al trabajar en dimensiones mas altas es conveniente utilizar el Transporte 
Paralelo de Levi-Civita (TP) para medir la curvatura. 

Supongamos dos puntos P y Q de M, y un arco a que los une. Entonces, 
el TP es un mapeo ortogonal del espacio vectorial S p en el punto P al espacio 
vectorial S q en el punto Q. En particular, si P = Q tenemos que a es un aro 
y el TP mapea S p en si mismo, por lo tanto TP debe ser una rotación de S p 
en un ángulo e(á) que depende del camino recorrido a. Este ángulo corresponde 
precisamente a la curvatura integral Gaussiana: 

e(a) = j K{P)dA (1.11) 

a 

ya que si por ejemplo la curva cerrada a corresponde a un trangulo t ten- 
dremos e(a) — EQ t . 

Además, e es una función aditiva en el camino que se recorre. Si tomamos 
dos aros a y b con el mismo punto final P y realizamos un TP alrededor del 
"producto" ab (i.e. recorremos primero el aro a y luego continuamos por el aro 
b), tenemos: 

e(ab) = s(a) + s{b) (1.12) 

No consideraremos casos en los cuales un vértice de M pertenece al aro a, 
ya que en este caso el TP no se puede definir sin ambigüedades. 

Si un aro o puede ser deformado continuamente hasta obtener un aro b, 
manteniendo P fijo y sin encontrarnos con un vértice, tenemos que e(a) = e(b) 
y por lo tanto escribiremos a = b. 



1.2.1. e— Conos 

Consideremos ahora una variedad M con un solo vértice. Esta variedad pue- 
de imaginarse entonces como un cono o una pirámide. Tomando coordenadas 
polares planas p y 6 con origen en el vértice, la métrica vendrá dada por: 

ds 2 = dp 2 + p 2 d6 2 (1.13) 

En el plano Euclideano dos puntos con coordenadas p iguales y con coorde- 
nadas 6 diferentes en un múltiplo de 2-k son idénticos. En la variedad M dos 
puntos con la misma coordenada p serán idénticos si sus coordenadas 6 difieren 
en un múltiplo de 2-k — e. 

Esta variedad es conocida como el e — Cono, y es Euclideana en todos sus 
puntos excepto en el értice, el cual tiene una deficiencia e. 
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Consideremos ahora la recta real R con variable z. El producto R®e — Cono 
es una variedad que posee la métrica: 

ds 2 = dz 2 + dp 2 + P 2 d0 2 (1.14) 

Este espacio es Euclideano excepto en la linea p = y es conocido como 
e — '¿Cono 

En general, si consideramos el producto R n ~ 2 ® e — Cono con la métrica 

ds 2 = dz\ + dz\ + ... + dzl_ 2 + dp 2 + p 2 d6 2 (1.15) 

obtenemos entonces un e — nCono y esta variedad es Euclideana excepto en 
el subconjunto plano n — 2 dimensional p — 0. 



Capítulo 2 

Variedades n-Dimensionales 



2.1. El Grupo Fundamental 

Consideraremos ahora una Variedad Rimanniana n-dimensional M. Tene- 
mos entonces las siguientes definiciones: 

Definicin 6 Una trayectoria a es una función continua y univaluada A(s) de 
un parámetro s € 1 , < s < 1. Los puntos A(0) y A(l) son los puntos inicial 
y final de la trayectoria. 

Si A(0) = A(í) = A la trayectoria se denomina Aro 

Definicin 7 Sean a y b dos trayectorias definidas por las funciones A(s) y B(s) 
tales que A(l) = B(0). Entonces el producto c = ab es definido por la función: 

r ¡ \ _ j A(2s) para s < 0,5 

Ü W-| B(2«-l) para s> 0,5 

Este producto de trayectorias es asociativo pero en general no es conmutativo 
(ab 7^ ba). 

Si a y b son aros con el mismo punto final P entonces ab y ba existen y son 
aros con punto final P. 

Definicin 8 Sea a una trayectoria definida por la función A(s). La trayectoria 
inversa a -1 esta definida por la función A(í — s). 

Consideraremos ahora que la variedad M es Euclideana excepto en un sub- 
conjunto cerrado w C M que posee curvatura. Sean entonces a y b dos trayecto- 
rias con los mismos puntos iniciales y finales, y definidas por las funciones A(s) 
y B(s). Ademas supondremos siempre que A(s) ^ w y B(s) ^ w. Definimos 
entonces: 

Definicin 9 Si existe una función univaluada P(s,t) des^R,0<s<ly 
deíeR,0<t<l continua en los parámetros t y s , tal que P(s,0) — A(s) 
y P(s, 1) = B(s) decimos que se ha deformado a en b y la función P(s, t) se 
llama una deformación. 
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Deflnicin 10 Sea una deformación P(s,t) que deforma a enb tal que P(s,t) £ 
w para todo s,t. Escribimos entonces a <~ b. 

Teorema 1 El símbolo ~ satisface las propiedades formales de una equivalencia 
(holonomia) . 

Para probar el anterior teorema sean o y b dos trayectorias. 

Teneemos que a ~ a ya que podemos tomar como deformación la función 



P{s,t) = A{s) 

Si a <~ 6 entonces P(s, 1 — i) deforma b en a y P(s, 1 — í) ^ io por lo tanto 



b <~ a 

Por último, si P(s, í) deforma a en 6 y Q(s, í) deforma 6 en c y ademas a <~ b 
, b ~ c, entonces la función : 



deforma o en c y ademas ñ(s, t) ^ w por lo que a <~ a 

Las anteriores definiciones se aplican de la misma forma a los aros con la 
condición adicional de que después de la transformación la nueva trayectoria 
sea un aro, esto es P(0, í) = P(l,t) para todo t. 

Deflnicin 11 Sean a y b dos aros y P(s,t) una deformación que lleva a en b. 
Si se cumple que P(0, t) = P(l, í) = P(0, 0) = P entonces escribimos a w b. 

El signo w se aplica únicamente a aros. Claramente a w b implica a ~ b 
pero la implicación inversa no es cierta. 

Teorema 2 El símbolo w define una equivalencia (homotopía en M — w) 

Para probar este teorema sean a y 6 dos aros y P(s, í) una deformación que 
llea a en b. 

Observemos que o w o ya que podemos tomar como deformación la función 



Si a rí 6 entonces P(s, 1 — í) deforma b en o y P(s, 1 — í) ^ w por lo tanto 



Por último, si a, b y c son aros y P(s, í) deforma a en 6 y Q(s, í) deforma 6 
en c y ademas a ~ b , b ~ c, entonces la función : 




para t < 0,5 
para t > 0,5 



P(s,t) = A(s) 




P(s, 2t) para í < 0,5 

Q{s,2t-1) para í> 0,5 



deforma a en c y ademas P(s, t) ^ w por lo que a w c. 
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Definicin 12 Defnimos ahora el aro unidad u en un punto P por medio de la 
función U(s) = P = const. 

Claramente para cualquier trayectoria a tenemos aa^ 1 = u 

Teorema 3 Sean a, a, b y b' aros tales que a w d y b w b' . Entonces ab w aV. 

Para la prueba del teorema sean P(s,t) la función que deforma o en a y 
Q(s,t) la función que deforma b en b'. Entonces la función 

„, P(2s,t) para s < 0,5 

[S ' ' ~ \ Q(2s-Í) para s > 0,5 

es una deformación de ab a aV , con ñ(0, t) = R(l, t) y además R(s, t) ^ w . 
Por lo tanto tenemos ab w ai)'. 

Si tomamos ahora un conjunto de aros con punto final P, podemos dividirlo 
en clases de equivalencia. Cada una de estas clases esta caracterizada por uno 
de sus elementos, y por lo tanto utilizaremos uno de sus aros para identificar la 
clase completa. 

Definimos entonces el GRUPO FUNDAMENTAL como el conjunto de todas 
las clases de equivalencia con el producto de trayectorias dado anteriormente. 
Para demostrar que se cumplen los axiomas de un Grupo veamos que: 

Si a y b son clases entonces todos los productos aV con a w a y b w b' son 
homotópicos y pertenecen a una nueva clase que definimos como el producto ab. 
Este producto de clases es asociativo. 

Si o es una clase y u es la clase unitaria entonces au = ua = a. 

A cada clase o asociamos una clase inversa a -1 tal que aa^ 1 = a _1 a = u. 

A cada aro a se le puede asociar una matriz ortogonal S(a) de la siguiente 
manera: 

Tomamos un vector A en el punto P y lo transportamos paralelamente 
(TP) a través de a .Existe entonces un mapeo lineal y que preserva la norma 
que relaciona las posiciones inicial y final del vector A. Si escogemos un sistema 
de referencia ortogonal local en P, este mapeo se puede representar mediante 
una matriz ortogonal S(a). 

Como caso particular, si nuestro sistema de referencia es localmente pseudo- 
Euclideano (e.g. Lorentziano), las matrices que represetan el TP son matrices 
de Lorentz. 

Ahora bien, dos aros homotópicos producirán la misma posición final del 
vector transpoprtado paralelamente, y por lo tanto serán representados por la 
misma matriz. Esto nos indica que la matriz ortogonal S(a) es función de las 
clases de Equivalencia, lo que implica que el conjunto de matrices ortogonales son 
una Representación del Grupo Fundamental en el Espacio Vectorial del punto 
P. 
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Para comprobar esto, veamos primero que a todos los aros v ~ u les corre- 
sonde la matriz unidad. Para ello sea T(s, t) la función de deformación que lleva 
v en u. El conjunto de todos los puntos T(s,t) es una superficie simplemente 
conexa S. Podemos tomar entonces en E y en una vecindad de S apropiada, un 
sistema de referencia globalmente Euclideano, y asi el resultado de un TP a lo 
largo de v calculado en este sistema de referencia es la identidad. 

Ahora bien, sean S(a), S(b),S(c) son las matrices correspondientes a los aros 
a, b, c y c — ab. Tomamos un vector V en P y lo transportamos a lo largo de 
a. Obtenemos entonces 5(a) V. Si transportamos ahora este vector a lo largo de 
b obtenenmos 5(6)5(a)V. Este resultado debe ser igual a transportar el vector 
original V a lo largo de c, es decir S(c)V = S(b)S(a)V. Ya que el vector V es 
arbitrario obtenemos S(c) = S(b)S(a). 

Por ultimo, si a w a' entonces podemos escribir a = va donde v w u. Tenemos 
entonces que S(a) = S(v)S(a), pero ya que S(v) — 1 obtenemos S(a) = S(a). 
Es decir que la matriz es función de la clase de equivalencia. 



2.1.1. Grupo Fundamental para los e — Conos 

Como un caso particular tomamos un aro en un e — Cono. Podemos entonces 
dar el ángulo 9 para cada punto en el aro como una funcioón de s. Podemos 
entonces definir sin pérdida de generalidad 9(0) = 0. De esta manera tenemos 
entonces 9(1) = N(2ir—e) donde N nos da el número de vueltas que el aro da a la 
linca p — (Hueso) el cual coincide con el subconjunto w definido anteriormente. 
Dos aros con igual N son homotópicos. De esta manera llamaremos a(N) la clase 
de equivalencia correspondiente. 

El producto de dos clases de equivalencia será a(N)a(M) — a(N + M). Es 
decir que el Grupo fundamental es Isomorfo al grupo de los enteros bajo adición. 

Para encontrar la representación del Grupo Fundamental en términos de 
las matrices Ortogonales S(a(n)) = S(n) notemos que si llamamos 5(0) = S 
tenemos que para un N dado S(N) = S N . 

S es llamdo el Generador del Hueso. 

Consideremos ahora un s — nCono. Tomamos un vector V en el punto P, y lo 
transportamos a lo largo de a(l). E lvector V puede expresarse en términos de 
dos componentes ortogonales, una en el e — Cono y otra en R n ~ 2 . La primera 
componente será rotada en un ángulo e, mientras que la segunda permanece 
invariante. Esto nos muestra que la representación del Grupo Fundamental en 
este caso es una matriz n x n con un espacio invariante R n ~ 2 . por ejemplo en el 
caso de un e — 3Cono el espacio invariante es una Recta y por lo tanto la matriz 
S es una rotación con eje R n ~ 2 = R. 
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2.2. Espacios Esqueleto 

Los Espacios Esqueleto corresponden a la generalización n-dimensional de 
los Poliedros definidos anteriormente. En ellos la curvatura se concentra en sub- 
espacios de dimensión n — 2 llamados Huesos. Una de estas estructuras son los 
e — nConos. 

Al igual que en el caso bidimensional comenzamos por una descomposición 
simplicial de la variedad, en este caso n-dimensional M. Esta descomposición 
determina la topología de la variedad, pero no su métrica, la cual será definida 
siguiendo los siguientes axiomas: 

Axioma 1 La métrica en el interior de cualquier simplex cerrado n-dimensional 
T n es Plana (e.g. Euclideana en M 3 o Lorentziana en M4). 

Podemos entonces calcular la distancia entre dos puntos cualquiera dentro 
del Simplex, definir un sistema de coordenadas cartesiano y expresar las coor- 
denadas de los puntos de la frontera de T n en este sistema. 

Axioma 2 En la métrica propia de T n , la frontera del Simplex se puede expre- 
sar como la suma de n + 1 simplex cerrados T„_i y estos simplex son planos. 

Axioma 3 Si el simplex T „_i es frontera común de T n y de T n entonces la 
distancia entre dos puntos cualquiera de T„_i es la misma en los sitemas de 
referencia de T n y de T n . 

Axioma 4 Si P € T n y V € T n y P,P" están lo suficientemente cerca de la 
frontera común T„_i entonces la distancia PP" se define como el ínfimo de las 
distancias PQ + QV para Q G T„_i. 

Con estos axiomas obtenemos una métrica plana también en T„_i excepto 
en las fronteras T n -2- Además esta definición es estrictamente local, pero nos 
permite unir suavemente los simplex vecinos. 

De esta manera, la Variedad M es Plana excepto en la unión w de los simplex 
frontera o Huesos T„_i. 

En particular en una variedad tridimensional M 3 los Huesos son segmentos 
de recta mientras que los T„_ 3 son puntos. Cada hueso conecta dos T n -3- Si 
conocemos la longitud de todos y cada uno de los huesos, entonces podemos 
conocer la estructura geométrica de los simplex y de la variedad. En particular 
podemos conocer también el ángulo dihedral. Por ejemplo, si en un hueso dado 
no existe curvatura, entonces la suma de todos los ángulos dihcdralcs alrededor 
del hueso debe ser 2tt. Cuando en un hueso existe curvatura, la suma de los 
ángulos dihcdralcs sera 2n — e. Reconocemos entonces que la sección transversal 
de este hueso sera localmente un e—Cono, y por lo tanto el hueso sera localmente 
un e — iCono. En general en una variedad n-dimensional M , el hueso T„_2 
tendrá localmente la geometría de un e — nCono. 
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Ahora bien, un T„_3 puede ser en general un extremo de m huesos T^_ 2 T^_ 2 T^ 
. Si estos ra se orientan de manera tal que tengan la dirección positiva saliendo 
de T„_ 3 entonces se dice que T„_ 3 es una m-juntura orientada. 

En una m-juntura orientada existe una relación entre los generadores de los 
huesos S p (p : 1, 2, m). Para encontrarla supongamos una juntura idealizada 
aislada ( i.e. los huesos se extienden hasta el infinito). Seleccionamos entonces 
un orden cíclico para los huesos. Sea A p el sector plano que tiene como lados los 
dos huesos contiguos T%_ 2 y T^_\. Supondremos que dos de estos sectores no 
tienen puntos en común, excepto tal vez los puntos de la frontera (hueso). De 
esta manera estos sectores dividen el espacio en dos regiones M'y M". Tomamos 
ahora un punto Pe M'y un punto Q G M'y m trayectorias t p que unen a P con 
Q. Supongamos ahora que la trayectoria t p intersecta al sector A p una y solo 
una vez. De esta manera t p t~_ 1 es un aro a p con punto final P y que encierra 
al hueso T%_ 2 una vez. Este aro es un candidato para ser el generador de T%_ 2 . 

El producto de todos los aros formados de esta manera es: 

010203. ...o p ¡=s u (2-1) 

o en términos de los generadores de los huesos: 

S V S 2 S 3 ...S P = 1 (2.2) 

Como veremos más adelante, la ecuación 2.1 corresponderá a la aproximación 
geométrica de la Identidad de Bianchi en una Variedad Diferenciable Riman- 
niana. 



2.3. Variedades Diferenciales 



Para realizar la transición del Espacio Esqueleto a una variedad Diferenciable 
incrementaremos la densidad p de Huesos manteniendo el producto pe como una 
función suave, tal como se describió anteriormente. 

Sea M3 una variedad tridimensional. Consideremos un haz de Huesos para- 
lelos en M3 . Supondremos entonces que la curvatura inducida por los huesos es 
pequeña y por lo tanto M3 es aproximadamente Euclideana.Sea U un vector 
unitario paralelo a los huesos. Para medir la curvatura de la variedad trans- 
portamos paralelamente un vector A alrededor de un pequeño aro de área S y 
normal unitaria n. Es decir que tenemos X = En. Después de realizar el TP el 
vector A ha rotado alrededor de U en un ángulo a = Ne donde N es el numero 
de Huesos acerrados por el aro. y esta dado por N = p(U ■ S). Asi, la diferencia 
entre el vector resultante al terminar el TP y el vector inicial será: 



5A = pe(U-E)(UAA) 



(2.3) 
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Ahora bien, de la Relatividad General tenemos que el Transporte Paralelo 
de un vector A^ produce el cambio: 

H« = Rl a[j ^A a (2.4a) 

donde S a/3 = e a/37 X¡ 7 con e al3j el simbolo de Levi-Civita en tres dimensiones 
y ^ es el tensor de Riemann. 

Comparando las Ecuaciones 2.3 y 2.4a obtenemos: 

Rficra/S = peUfitjUafi (2.5) 

Con U^a = e ¡llJ \U x . Esta aproximación tiene todas las propiedades de si- 
metría que posee el Tensor de Riemann, y ademas cuando existen más haces de 
huesos, el efecto es puramente aditivo. 

En una variedad n-dimensional M, un hueso tiene una orientación detere- 
minada por un tensor Sesqui-Simétrico U^ a tal que U^U^ = 2 y U^crU a p + 
U^aUfju + UfjtfiUua = 0. En este caso, si trabajamos con un haz de huesos parale- 
los, escogemos las coordenadas x\ y X2 perpendiculares al haz y las coordenadas 
£3,2:4, ...,x n paralelas al haz. Así, tenemos U12 — —U21 — 1 mientras que las 
demás componentes son cero. La Ecuación 2.5 aún es válida si hacemos que los 
Índices /1, v, a, (3 tomen los valores 1, 2, 3, n. 

Por otro lado, independientemente del valor de n, el escalar de curvatura 
viene dado por: 

R = 2pe (2.6) 



2.3.1. Identidades de Bianchi 

Mostraremos ahora como las Ecuaciones 2.1 y 2.2 corresponden a las Iden- 
tidades de Bianchi. Para ello observemos que la Ecuación 2.1 nos indica que al 
recorrer todos los aros alrededor de los huesos en una m-juntura en M 3 obtene- 
mos el mismo vector con que comenzamos, es decir que la composición de los 
m aros produce la identidad. De esta manera, si las deficiencias de cada hueso 
de la juntura e son pequeños, podemos escribir esta relación para la m-juntura 
como: 

m 

E £ P ^ = ( 2 - 7 ) 

p=l 

o equivalentemente 



ra 

^2s p U p = 



(2.8) 
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Ahora bien, supongamos una distribución de densidad g de m-junturas idén- 
ticas. Es decir que tendremos m haces de huesos paralelos. Supongamos entonces 
que el p-esimo haz de huesos tiene una deficiencia e p , una dirección dada por 
el vector U p y una densidad p p . Es claro que p p no es constante, ya que cada 
juntura hace un decaimiento o un crecimiento de esta densidad con una razón 
que depende de la densidad de junturas g. 

Por definición tenemos que p p es la densidad de flujo de huesos que salen 
de la juntura a través de una superficie X ortogonal a U p . Supongamos ahora 
que la posición de la superficie S esta determinada por una abscisa s que se 
mide sobre la dirección del vector U p . Sea entonces C un cilindro de alura ds 
y base S(s). Así, cada juntura dentro de C incrementa en una unidad el flujo 
de p p huesos a través de la superficie E(s + ds). Además dentro del cilindro C 
existirán £ • g ■ ds junturas. De esta manera tendremos: 

Ep p (s + ds) - Spp(s) = llg(s)ds (2.9) 

es decir: 

En general, esta expresión será: 

U*'" ■ 0„ft, = <? (2.11) 

Ahora bien, el Tensor de Riemann vendrá dado por la expresión 2.5 pero 
debemos sumar sobre todas las junturas, es decir: 

R af 3X8 = Y. £ PPP U lp U ls ( 2 - 12 ) 



Ya que hemos supuesto que la curvatura inducida por los huesos es pequeña, 
podemos aproximar también las derivadas covariantes con derivadas normales, 
y de esta manera las identidades de Bianchi serán: 



d d d 

B af ¡XS 7 = Q-^R a f3XS + -g^jRapSy + g^gR a /3 7 X = (2.13) 



Reemplazando la Ecuación 2.12 en la Ecuación 2.13 y tratando e p y U p como 
constantes tenemos: 
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dx x 



"dx 5 



es decir: 



B a/ 3\s- ( - ^SpU^p 



JL + U P JL + U r> A- 

dxt & ~<dx x ~< x dx s 



Pp = O (2.15) 



Es inmediato verificar que se tiene la identidad entre operadores: 



d d d d 



d d d 

U »to +U *'dxX +U *d¡? =e ->» ÜVa ° 



(2.16) 
(2.17) 



Reemplazando esta identidad en 2.15 tenemos: 



B al 3\8y = ^SpUlp [e 1 x 5 U p ' v d v \ p p = 
p 

Bapxsy = e 1 x¡^£ p Ul [} U p ^ 'd v p p = 



(2.18) 
(2.19) 



Haciendo uso de 2.11 tenemos: 



(2.20) 



que coincide exactamente con 2.7 y por lo tanto nos muestra la relación de 
2.1 y de 2.2 con las identidadesde Bianchi. 
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2.4. Variedades Lorentizianas 

Cuando la signatura del tensor métrico g de una variedad n-dimensional 
M n , es sig(g) = 2 — n la variedad se denomina Lorentziana. En este caso, la 
métrica de la variedad no es siempre positiva. En este caso, la representación del 
grupo fundamental, es decir las matrices ortogonales 5* deben ser reemplazadas 
por matrices de Lorentz. Correspondientemente, en el Espacio Esqueleto, los 
ángulos de déficit pueden ser imaginarios. 

Nos restringiremos ahora al caso 4-dimcnsional. En esta variedad, los huesos 
son triángulos, y la métrica tiene una coordenada temporal. De esta manera 
tenemos tres clases de huesos: 

1. Como- de-Espacio: Cuando cualquier vector en el hueso es como de espa- 
cio. El ángulo de déficit e es imaginario. 

2. Nulo: Cada vector en el hueso es una combinación lineal de un vector nulo 
y un vector como de espacio ortogonal a él. En este caso tenemos s = . 

3. Como-de-Tiempo: Existen vectores como-de-tiempo y como-de-espacio en 
el hueso, e es real. 

Ahora definiremos el área L de un hueso. Para ello, notemos que en M 3 el 
área de un triángulo que tiene como lados los vectores B y C viene descrita por 
el vector: 



2A = BxC (2.21) 
que es perpendicular a B y a C. La magnitud del área sera: 

AA 2 = B 2 C 2 sin 2 8 (2.22) 
donde 8 es el ángulo formado por los vectores B y C. Tenemos entonces: 

AA 2 = B 2 C 2 (l - eos 2 8) = B 2 C 2 - B 2 C 2 eos 2 8 (2.23) 

AA 2 = B 2 C 2 - (B • C) 2 (2.24) 

En el espacio 4-dimcnsional, definiremos entonces el área L como la genera- 
lización: 
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4L 2 = (B^ f - (BUBJiCCn) (2.25) 

Se puede comprobar entonces que al igual que en el caso 3-dimensional, el 
vector 2L es dual al bivector construido de B a C. Es decir que si B va en 
dirección i° y C ra en dirección a; 3 , entonces L es un bivector en el plano 
(x 1 , x 2 ). Asi, L es una cantidad real siempre que B y C no sean ambos como- 
de-espacio. Por lo tanto trabajaremos únicamente con huesos nulos y como-de- 
tiempo. Además es claro que si el hueso es nulo, tenemos L = 0. 
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Capítulo 3 



Cálculo de Regge y 
Relatividad General 



3.1. Ecuaciones de Campo de Einstein en el Va- 
cio 

Utilizando todas las anteriores definiciones y el esquema anterior del Cálculo 
en Espacios Esqueleto, deduciremos ahora el análogo discretizado de las Ecua- 
ciones de Campo de Einstein para el Vacio. 

La manera más directa de llegar a la esqueletización de un espacio de Eintein 
y de llegar a las ecuaciones de campo es partiendo de un principio variacional. 
Para ello tomamos la Acción de Einstein- Hilbert: 



Las ecuaciones de de Campo de Einstein corresponderán entonces a la con- 
dición de que / sea un extremo. Es decir 51 = 0. Esta condición es aplicable 
solamente cuando consideramos un espacio libre de materia y de campo electro- 
magnético, es decir que obtendremos las ecuaciones de campo para el espacio 
vacio. 

En el espacio esqueleto, la variación de esta acción debe realizarse sobre las 
longitudes de los simplex Ti, ya que como sabemos estas longitudes contienen 
la misma información que el tensor métrico. 

Ahora bien, para obtener una expresión de / en términos de las longitudes, 
observemos primero que ni T4 ni T 3 contribuyen a la integral, ya que R es sim- 
plemente una función distribución con soporte w tal como vimos anteriormente. 




(3.1) 
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Por lo tanto / debe ser una función de los huesos T^. Si k es el número de 
huesos en el esqueleto tenemos I = J2 F ( T 2) ■ Ahora bien, el área del hueso 

k 

L es una función aditiva, mientras que la función F es la misma para todos 
los huesos. Supondremos entonces que F es proporcional al área del hueso, i.e. 
F(T|) = L k f(e k ). 

Por otro lado, un hueso T% puede ser expresado como la superposición de 
dos huesos T 2 fe 'y T 2 ' £ ' con I a misma forma y área y tales que e k = e k + £k ■ Se 
sigue entonces que /(efe) = /( £ fc) + /( £ *0- Esta condición implica que f(e) = Ce 
donde C es una constante. Así tenemos que la acción se puede expresar como: 



/ = Cj2L k £k (3.2) 



Utilizando la expresión para el escalar de curvatura dado por la ecuación 2.6 
y en los límites de una variedad aproximadamente plana tenemos que la acción 
viene dada por: 



1=^-1 egd 4 x (3.3) 



8tt 



Comparando las ecuaciones 3.2 y 3.3 tenemos C — g^, con lo que obtenemos 
finalmente la Acción de Regge: 



1 = h^ Lk£k (3 - 4) 

fe 



En una variedad n-dimensional M n tendremos la generalización: 

donde LJJ es la (n — 2)— dimensional medida del hueso T„-2- 

La siguiente tarea será encontrar la variación de la Acción de Regge en 
términos de las longitudes l p de los simplex Tf constituyentes de los huesos. 



Esto es: 
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SI = ^SY^Lkek = O (3.6) 



La dependencia de el área del hueso en esta longitud es clara, pero la depen- 
dencia del ángulo de déficit e es bastante complicada. Sin embargo, siguiendo a 
T. Regge'l' se puede comprobar (ver Apéndice 1) que la variación puede realizar- 
se como si los e n fuesen constantes, es decir que la ecuación 3.6 es simplemente: 




En el caso de M 3 los huesos serán los vértices y por lo tanto SL k = Slk, y 
por lo tanto la Ecuación 3.7 se reduce simplemente a la condición = 0. Esto 
quiere decir que todos los espacios esqueletos de Einstein de dimensión 3 son 
planos. 

Ahora bien, en una variedad 4-dimcnsional, y de acuerdo con la trigonometría 
elemental, el cambio en el área de un triángulo en función de la variación en la 
longitud de sus lados será (ver Apéndice 2): 

SL k = ^y^¿ p cot (6p k ) Sl p (3.8) 
v 

Donde 6 p k es el ángulo opuesto al lado l p en el triángulo T 2 fc . De esta manera 
reemplazando este resultado en la Ecuaciñon 3.7 obtenemos: 

si = -^J2J2 £kl p cot ^ si p = ( 3 - 9 ) 

fe p 

es decir: 

^£fc cot (9 p k) = con p=l,2, .... (3.10) 

fe 

La ecuación 3.10 corresponde a la versión discretizada de las Ecuaciones de 
Campo de Einstein en el vacio según el método de Regge. Se puede observar 
que se tiene una ecuación por cada longitud de un vértice dentro de la región 
de espacio-tiempo que se analiza. 
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3.2. Desviación de la Geodésicas en el Modelo 
de Regge 

Una de las preguntas abiertas en la formulación de Regge es la interacción 
entre el campo Gravitacional y materia. Se han logrado realizar la inclusión de 
algunos términos que representan materia en esta formulación, pero para ello se 
deben hacer algunas restricciones acerca de la clase de materia que se acopla al 
campo gravitacional y se deben imponer ciertas simetrías al campo. 

Otra manera de incorporar la materia en la formulación de Regge es la 
de identificar primero los parámetros que identifican a la materia tanto en la 
formulación continua de la Relatividad General como en el cálculo de Regge, 
y de esta manera simplemente se insertarían los términos apropiados en las 
ecuaciones discretas. Sin embargo esta identificación no es trivial y usualmcntc 
conduce a divergencias. 

Una manera de comenzar con esta identificación es analizar inicialmente 
una de las manifestaciones mas elementales de la influencia de la gravedad en la 
materia: la Desviación de Geodésicas. Esta desviación nos describe la aceleración 
relativa de partículas de prueba libres en un espacio-tiempo curvo. Primero 
analizaremos el caso en la formulación continua siguiendo a Synge^ . 



3.2.1. Desviación de Geodésicas en Relatividad General 

Sea C[v) un conjunto de curvas con ecuación — x^(u,v) , de tal manera 
que v es constante a lo largo de cada curva, y parametrizada por u. Este conjunto 
de curvas forman una superficie 2-dimensional. 




Figura 3.1: Familia de Geodésicas C(v) parametrizadas por u. Los vectores U 
y V son tangentes a las curvas con v y u constante respectivamente. 
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Los vectores = d u x» = y — dyX» — satisfacen: 

SU» SV» 



Sv 

S 8 



(3.11) 



donde j¿ y j¿ son las derivadas covariantes a lo largo de las curvas con v y 
u constante respectivamente. 

Consideremos dos curvas adyacentes C(v) y C(v + Sv). De esta manera, el 
vector V* representara la separación entre las dos curvas. Para encontrar como 
se desvia la curva C(v + Sv) de la curva C(v) calculamos la segunda derivada 
covariante del vector de separación, a lo largo de la curva C(v), parametrizada 
por u. Esto es: 



s 2 v» s sv» (3 12) 



Su 2 Su Su 



Utilizando 3.11 tenemos: 



s 2 v» = ±s_u»_ 

Su 2 Su Sv 

Ahora bien, sabemos que la diferenciación covariante en una 2-supcrficic 
satisface la regla de conmutación: 



Así tenemos que la ecuación 3.13 será: 

^ = ^---^ 

donde R^p^ es el tensor de Riemann. 

Si tomamos las curvas C(v) como geodésicas tenemos la condición: 
Así, tenemos que la ecuación de la desviación de las geodésicas será: 



5 2 V» 
Su 2 



S 2 V» 
Su 2 

K Pl u a u vi 



Kfh uauf>vr 



(3.17) 
(3.18) 
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Además, utilizando las simetrías del Tensor de Riemann tenemos: 
5 2 V^ 



-j^^K^v 13 ^ = ° ( 3 - 19 ) 

Ahora bien, si en lugar del parámetro u utilizamos la longitud de arco s, y 
las lineas parametrizadas con v las tomamos ortogonales a las C(v); la ecuación 
de desviación de Geodésicas toma la forma general: 



^ + R(U,V,U)=0 (3.20) 

Como caso especial, en una variedad 2-dimensional el tensor de Riemann 
tendrá únicamente una componente diferente de cero, que tomaremos como 



3.2.2. Desviación de las Geodésicas en el modelo discreti- 
zado: 



Para observar el efecto de la desviación de las geodésicas en el caso discreto 
seguiremos el trabajo de S. Chakrabarti, et. al.^1. En este caso será suficien- 
te considerar un espacio-tiempo 2-dimensional. De esta manera la curvatura 
estará concentrada en los vértices de la red simplicial. 




(a) (b) 



Figura 3.2: (a) Dos geodésicas originalmente paralelas, convergen cuando entran 
al triangulo [023] debido al ángulo de déficit e en el vértice [0] . (b) Al realizar un 
corte en la arista [02] y colocando el simplex en un plano se observa el ángulo 
de déficit e. Las geodésicas son ahora rectas, y no cambian de dirección al pasar 
de un triangulo al otro. La desviación solo se hace visible en el triangulo [023]. 



En la figura 3. 2 (a) vemos como dos geodésicas que son paralelas al entrar en 
el triángulo [014] sufren una desviación debido al ángulo de déficit e del vértice 
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[0]. Para obtener una mejor visión de lo que sucede en realidad realizamos un 
corte en una de las aristas del 3— simplcx [01234] y colocamos la figura resultante 
en un plano (Figura 3.2(b)). Así, las geodésicas se convierten en simples rectas. 
Sin embargo, aún cuando inicialmcntc son paralelas, al llegar al triángulo [023] 
convergen debido al ángulo de déficit e. 

Observamos entonces que la pasar alrededor de un vértice, las geodésicas 
sufren una desviación que se manifiesta en una cambio en el ángulo de con- 
vergencia de las geodésicas Aa , y este cambio corresponde al ángulo de déficit 
asociado con el vértice, i.e. 



Es decir que al igual que en el caso continuo, la curvatura enfoca o dispersa 
las geodésicas. Sin embargo, en el caso de un único vértice, el cambio en el 
ángulo de convergencia de las geodésicas no puede ser continuo. Así, no es 
posible encontrar una ecuación de desviación de geodésicas que dependa de la 
segunda derivada del vector de separación, tal como en el caso continuo. 

Desviación Producida por una Distribución de Vértices (Caso 2-dimensional) 

Para lograr encontrar el análogo de la Ecuación de Desviación de Geodésicas 
(Ec. 3.19 ) en la formulación de Regge, estudiaremos el caso de un espacio 2- 
dimcnsional. 

Consideramos una región de una red Simplicial, en la cual existen N vértices 
con ángulos de déficit {si\f =1 ■ Asumimos además que cada ángulo de déficit e¿ 
es pequeño, y que el área que contiene los vértices, AA, es también pequeña. 
Sin embargo, el ángulo total de déficit por unidad de área 



es finito. Además hemos definido n como la densidad de vértices por unidad 
de área y e como el ángulo de déficit promedio por vértice: 



Aa = e 



(3.21) 




(3.22) 




(3.23) 



Consideremos ahora dos geodésicas C\ y C*2 que pasan rodeando el conjunto 
de N vértices. El cambio en el ángulo de convergencia de las dos geodésicas es- 
tará determinado por el ángulo total de déficit asociado con el área considerada: 
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N 

Aa = (3.24) 

¿=i 

Asumiremos ahora que cada geodésica esta parametrizada por la longitud 
de arco s y que estas geodésicas están definidas además por un parámetro v tal 
que las lineas s = const. son ortogonales a las geodésicas. Aproximamos ahora 
esta situación a una variedad 2-dimensional con coordenadas (x 1 ,^ 2 ), tales que 
las lineas con coordenada a: 1 = const. corresponden a las curvas con v = const. 
mientras que las lineas con x 2 = const. corresponden a las geodésicas. 




Figura 3.3: (a) Las geodésicas C\ y C2 son defelctadas por N vértices, (b) 
Representación de la desviación de geodésicas producidas por la distribución de 
N vértices. 



De esta manera, el área que contiene los N vértices esta dada por AA = 
x 1 Ax 2 m. (Figura 3.3 (b)). 

Las componentes del Vector de Separación V y del vector velocidad de la 
geodésica U serán: 



V = (x\0) (3.25) 



i.í- 



U=(0,— ) =(0,1) (3.26) 



Además, de la Figura 3.3(b) podemos observar que el ángulo de convergencia 
se puede aproximar por: 



Ax 1 

a w tana = — — ¡r (3.27) 
Ax ¿ 
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De esta manera, utilizando las expresiones 3.27 y 3.22 podemos escribir la 
ecuación 3.24 como: 



Ax^ 



A (2?) = 5> = nmA (3 - 28) 
A (tí) = (3 ' 29) 



Ya que Ax 2 — As y como la norma del vector de separación es V = x 1 
tenemos: 



A (^j) + neVAs = (3.30) 

á(^r) +n ^° (3 ' 31) 



Llevando esta expresión al caso infinitesimal, y ya que el vector separación 
puede expresarse como V = (V^O) obtenemos: 



5 2 V X 

-jj- + neV 1 = (3.32) 

Ahora bien, como dijimos en la sección anterior, en el caso 2-dimcnsional el 
tensor de Riemann tiene una sola componente diferente de cero, que tomaremos 
como R\i2 , y de esta manera la Ecuación 3.20 puede escribirse como: 



S 2 V 1 

-—+^2^ = (3.33) 
os z 

Comparando estas dos ultimas ecuaciones obtenemos la identificación: 



#212 = ne (3.34) 

es decir que la única componente del Tensor de Riemann en el caso 2- 
dimensional se puede interprestar como el ángulo de déficit promedio por unidad 
de área. 
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Generalización a 4 dimensiones 

Para generalizar los anteriores resultados a una variedad 4-dimensional, sim- 
plemente debemos notar que en este caso la curvatura esta concentrada en 
2— simplex, es decir triángulos. Además, para una variedad 4-dimensional el 
tensor de Riemann tiene 20 componentes independientes. 

Si tenemos una región infinitcsiomal con N triángulos, una geodésica que 
se encuentre dentro de un plano ortogonal a uno de esos triángulos sufrirá una 
desviación igual a la magnitud del ángulo de déficit correspondiente al triángulo. 
Es decir exactamente igual al caso 2-dimensional. 

Ahora bien, en general la geodésica no va a encontrarse en el plano ortogonal 
a todos los triángulos, y además no todos los triángulos tienen la misma orien- 
tación. Si definimos el ángulo 9 como el formado entre la dirección del vector 
tangente a la geodésica y el vector normal a uno de los triángulos, tenemos que 
el vector tangente sufrirá una rotación en un ángulo e eos 9. 

De esta manera es inmediato que para hacer la generalización a 4 dimensiones 
solamente debemos reemplazar n por el número de triángulos por unidad de 
volumen y cambiar el ángulo déficit promedio e por el ángulo de déficit efectivo 
que afecta una geodésica con vector tangente U definido como: 



donde 0, es el ángulo entre el vector tangente U y el plano ortogonal al 
triángulo con ángulo de déficit e¿. 

De esta manera recuperamos la definición del tensor de Riemann obtenida 
en la Ecuación 2.5, es decir que en la formulación de Regge esta representado 
por el ángulo de déficit promedio efectivo por unidad de volumen. 



N 




(3.35) 



¿=i 



Apéndice 1: Variación de la 
Acción de Regge 



En este Apéndice mostraremos como para una variedad ra— dimensional M n 
, se cumple que para variaciones infinitesimales de l : 

J2 l p 6£ p = ( 36 ) 
p 

donde L™ es la medida (ra — 2)— dimensional del p— esimo hueso. 

Sea T n un simplex, y definamos en el un conjunto de coordenadas cartesianas. 
T n posee n + 1 simplex frontera denotados T£_ 1 (r = l,2,...,n+l). Dos de estos 
simplex T^_ x y T*_ x tendrán una frontera común denotada por T£f_ 2 - 

Sea XJ T el vector unitario normal al T^_ 1 (r = 1, 2, ra + 1), es decir 
(U r )^ (U r Y = 1. Si es el ángulo entre T^ l _ 1 y T*_ x entonces tenemos: 

C os6rs = {U r ) fl {U s Y (37) 

Ahora consideremos el tensor Sesqui- Simétrico {U r3 ) tiv = — {U rs ) vfí defi- 
nido por: 



sin (0„) (IT-)^ = (ET)^ - (ET)„ ([/«)„ (38) 

Claramente tenemos que (U 7 " 3 )^ (JJ rs ') tíV = 2 , por lo que este tensor es el 
mismo que se utilizó en la Sección 3.7 para calcular el Tensor de Riemann. Este 
vector define la orientación de T™_ 2 . 

Ahora bien, si L rs es la medida del simplex T^f_ 2 tenemos la identidad: 
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Y,L rs U rs = O 



(39) 



Para comprobar 39 tomemos por ejemplo un sistema en el que (U r ) — 
1 y {U r ) , Q = 0. Entonces las únicas componentes no nulas de (U rs ) vf ¿ son 
las que tengan v Y 0, y por lo tanto llamamos (V s ) v = (U rs ) v0 . Tenemos 
entonces (V s ) v (V s ) v = 1. Ahora bien, en T^_ l5 la cantidad £Z, rs (V s ) v A v es 

s 

simplemente el flujo del vector constante A v a través de la frontera de T£_i y 
por lo tanto debe ser cero. Así llegamos a la ecuación 39. 

Diferenciando la ecuación 37 tenemos: 



50 r 



sin 6 r 



(un,s(ur + (u s ),Hu r Y 



(40) 



Por otro lado, multiplicando 38 por (U r ) v obtenemos: 



<u rs ) (u r ) v = — - 

K Jv » y ' sin(9 r 

(u rs )(ur- 1 



(u r ),(u s ) v (u r r-(u r )Au r T(u s ), (4i) 



sin 6 r 



cos9 rs (U r )^-(U s ), 



(42) 



Multiplicando esta última expresión por S {U r Y y utilizando {U r )^ 5 (U r Y = 
tenemos: 



blli \Jf - 



(43) 



Reemplazando esta ecuación en 40 tenemos: 



se rs = (u rs ) vu WY s (u r Y + (u s Y s (u s Y] (44) 



Multiplicando por L rs y sumando obtenemos: 



J2L rs se rs = Y^Lrs {u rs ) vlí WY s (u r Y + (u s Y s (u s Y\ (45) 
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Utilizando 39 se tiene: 

J2 L rsSe rs = (46) 

r.s 

Si sumamos esta expresión sobre todos los simplex del esqueleto obtenemos 
finalmente: 

J2 l p S£ p = ( 4? ) 

V 

Es decir que la variación de la Acción de Regge sera simplemente: 

SI = -^S ^e k L k = g^X! \- LkÓ£k + £ kSL k ] = -^^2 £ kSL k = (48) 

fe fe fe 
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Apéndice 2: Cambio en el 
área de un triángulo en 
función de la variación en la 
longitud de sus lados. 

Consideremos el triangulo ABC de la Figura 4. Si hacemos una pequeña 
variación 6 A en uno de sus lados manteniendo constantes los lados B y C ob- 
tendremos el triangulo ABC que se observa. 




Figura 4: Cambio en el área de un triangulo 

El área del triángulo ABC será simplemente \Bh. Por otro lado el área del 
triángulo ABC sera \BK. De esta manera el cambio en el área es: 



SL= -B(h'-h) (49) 

De la Figura 4 observamos ademas que h — C sin 9 mientras que K = 
Cs'mO'. Asi tenemos: 
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8L=^BC{smff-sm6) (50) 

Utilizando el Teorema del Coseno tenemos A 2 = B 2 + C 2 - 2BC eos 9. De 
donde, diferenciando, obtenemos: 



2ASA = 2BCsin{9)89 (51) 
ASA 



BCsin9 



66 (52) 



Ahora bien, el ángulo #'será simplemente ff = 9+89. De esta manera tenemos 
que: 

sin ff = sin(6> + 89) = sin 9 eos 89 + eos 9 sin 89 (53) 

Ya que estamos considerando 8 A infinitesimal, tendremos también 89 infini- 
tesimal y por lo tanto podemos aproximar: 

sin ff = sin 9 eos 89 + eos 9 sin 89 « sin 9 + cos(9)S9 (54) 
Reemplazando en la ecuación 50 tenemos: 

SL = ^BC (sin 9 + cos(9)S9 - sin 9) = ^BC cos(9)S9 (55) 

Utilizando ahora la ecuación 52 obtenemos: 

1 / A8A \ 

8L = -BC eos 9 — - - (56) 



2 \BCsin9 
SL = ^ cot(9)A8A (57) 

Si ahora numeramos los lados del triangulo como li,h,h y realizamos va- 
riaciones a cada uno de ellos tendremos que la variación en el área del triángulo 
será: 



SL = ^J2 l P cot ( e p) Sl P ( 58 ) 
p 

donde 9 P es el ángulo opuesto al lado l p . 
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